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Abstract 

A new, more efficient approach to include a three nucleon force into three- 
nucleon continuum calculations is presented. Results obtained in the new and 
our old approach are compared both for elastic nucleon-deuteron scattering 
as well as for the breakup process. The advantages of the new scheme are 
discussed. 
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Three-nucleon forces (3NF) are getting more and more into the focus of few-nucleon 
studies, both theoretically [[I]] and experimentally ||. Therefore optimal forms of equations 
should be used to handle those forces most efficiently 

One of the first 3N continuum calculations including a 3N force were performed in || . 
Separable nucleon-nucleon (NN) forces together with a ir — ir exchange 3NF [Q averaged 
over spin and isospin degrees of freedom have been used. Thereby the standard Lovelace 
equations || have been generalized to include a 3NF. Starting from the AGS form fTj of the 
arrangement operators for identical particles and including a 3N force one gets [§] 

U = PG 1 + (1 + P)t A + PtG U + (1 + P)UG tG U (1) 

Here Go is the free 3N propagator, t the NN t-matrix, t 4 the corresponding quantity now 
driven by the 3NF V4 alone, and P the sum of a cyclical and anticyclical permutation of three 
objects. That operator relation (H) has to be applied on a channel state \<p > corresponding 
to Nd scattering, which is a product of the deuteron wave function and the momentum 
eigenstate of the projectile. Once U has been determined one gets the operator for the 
breakup process by quadrature. For more details and our notation in general we refer to ||. 

Eq. (|1|) has been solved in |3|] assuming NN t-matrices of finite rank. For the sake of a 
simple notation we just take the schematic form 

t=\g>r<g\ (2) 

which converts ([I]) into the structure of a single particle Lipmann-Schwinger equation 

X = [Z (2) + Z (3) ) + [Z (2) + Z^)tX (3) 

where 

X = < g\G UG \g > (4) 
= < g\PG \g > (5) 

and 
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Z 0) = 1 < ^| Go(1 + p)i 4(1 + P)Go |^ > (6) 

In addition to the well known particle exchange term there is now a potential term 
driven by the 3NF. We used 

(1 + P)U = U(l + P) = ~(1 + P)U{1 + P) (7) 

In ref. || Eq.@ has been solved under the simplifying assumption t 4 = V4. We would like 
to mention another approach in the context of finite rank NN forces 0]. There a simple ad 
hoc ansatz for £4 has been chosen. 

For general forces, which are not of finite rank, another set of equations within the 



Faddeev scheme was formulated |8|-[T0| and has been heavily used in |TT|-[T3|. This is a set 



of two coupled equations in case of identical particles. They have the form 

T = tP + tG T 4 + tPG T (8) 
T 4 = (1 + P)U + (1 + P)t 4 G T 

Once T and T 4 are found the operator for the 3N breakup process is given as 

U = (1 + P)T + T 4 (9) 

and the operator for elastic Nd scattering as 

U = PGq 1 + PT + T 4 (10) 

All these operator relations have to be applied onto the initial channel state \<p>. 

We have solved that coupled set |T]-|13|] by expanding t 4 , T and T 4 in powers of V4, which 
leads to a recursive set of equations . This is feasible and numerically precise but requires 
extensive resources in storage and computer time. 

For the 3N bound state in a momentum space treatment an algorithm has been proposed 
14jl , which has been taken up again in |y|. There the property of 3NF's applied up to now 



have been used, that it can be split naturally into 3 parts, each one of which is symmetrical 



under exchange of 2 particles, like the pair forces. For the ir — ir exchange 3NF for instance 
that splitting occurs automatically, since there are three possibilities for choosing a nucleon 
which undergoes the (off-shell) n — N scattering. Using that decomposition one can combine 
each of the three NN forces with a corresponding part of the 3NF having the same symmetry 
under particle exchanges. Then the derivation of the Faddeev equation from the Schrodinger 
equation follows exactly the same line as for NN forces alone and one arrives at one Faddeev 
equation using the identity of the particles. For the 3N bound state it has the form 

if> = G Q tP^j + G (l + tG )vi 1] (l + P)if> (11) 



This equation has been used in [|Tf| and later work, for instance in [[16 . 

Based on that formal insight it is obvious that a corresponding form should exist for 3N 
scattering. The elaboration of that expectation is the content of the present work. 

Let us split the 3NF into 3 parts 
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U 4 = E^4 (12) 

i=l 

as described above and let us arrange the 3N Schrodinger equation in integral form as 



y = G j2(V> + v} l) )y (13) 
i=i 

We used the usual "odd man out" notation for the NN forces. This is a correct homogeneous 
integral equation for a scattering state, which is initiated in the nucleon-deuteron (Nd) 
channel ||. Then splitting \l/ into three Faddeev components ipi one S e ^ s f° r the first 
component, for instance, 

V>i = GoiVt + Vf } )(^ + ^ 2 + V> 3 ) (14) 

The NN t-matrix t\ related to V\ is then introduced in the normal manner and one ends up 
easily with 

ii>i = <k + Goh(?p2 + fa) + (1 + G ti)G U 4 (1) (^i + ^2 + $j) (15) 
Here 0i is the initial channel state, as described above. 



For identical particles 

^ 2 + ^ 3 = P^i (16) 

and dropping the index 1 we get 

if) = + GotP^j + (1 + G Q t)G vi 1] {l + P)ij (17) 

This is the Faddeev equation for 3N scattering initiated in a Nd channel. 

Using Got = GV, where G = (E — H — V)^ 1 , one can easily extract the asymptotic 
behaviour in the elastic and breakup channels ||17|| . Thus the amplitude accompanying the 
familiar outgoing wave in the hyperradius is 

f = (1 + tG )VPij + (1 + tG )V 4 (1) (l + P)if> 

= tPij+(l + tG )v} 1 \l + P)ij (18) 

The full breakup amplitude constructed from all three Faddeev amplitudes is then 

U = (l + P)f (19) 

Now we see directly comparing and (|18D that 

ifj = 4> + G Q f (20) 

and therefore 

f = tPcj) + (1 + tG )vf } (1 + P)4> + tPG f + (1 + tG )vf } (1 + P)G T (21) 

This is the searched for single equation, which replaces the coupled set @. For = it 
reduces to the form, which we always use in solving the 3N continuum for NN forces only 
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The asymptotic behaviour of fllTD in the elastic channel yields the operator for the elastic 
process 

U = VPip + V} 1] (1 + P)ifj (22) 



The elastic amplitude arises by projecting U from the left with the channel state < <f>\, 
which carries the correct energy. Then V in the driving term can be replaced by Gq 1 . Using 
that and the relation ( pop we find 

U = PG 1 + PT + V} 1] (1 + P)4> + V A {1 \l + P)G f (23) 

This is the equation which provides U once T has been found. 

It remains to exhibit the formal connection to the previously used forms (|8|)-(|T0|). 

We use now the following property of V4 underlying its decomposition used throughout: 

(1 + P)V 4 = (1 + P)F 4 (1) (1 + P) (24) 

Let us then define 

V4 = U 4 (1) (1 + P) (25) 



and 



% + Wok (26) 



Then it follows 



1 + P)U = (1 + P)[V A + V 4 G V 4 + V 4 G V 4 G V 4 + ■■■} 

= (1 + P)V A {l \l + P) + (1 + P)K 4 (1) (1 + P)G V 4 + ■■■ 
= (1 + P)[V4 (1) (1 + P) + V} 1] (1 + P)G V 4 {1) (1 + P) + • • ■ 
= (l + P)[V 4 + V 4 G'oV r 4 + ---] 

= (1 + P)U (27) 



Further we introduce 



r 4 = t 4 + UGqT 

= V 4 + V 4 G n + V 4 G T (28) 



Applying (1 + P) from the left we get 



(1 + P)r 4 = (1 + P)U + (1 + P)UG Q T (29) 
= T 4 

according to (§). 

Let us now add the information from the first equation of the set (|8|) to get 

T + n = tP + tG (l + P)r 4 + tPG T + V 4 + V 4 G t 4 + V 4 G T (30) 

Using again (|28|) yields 

f = tP + V 4 + tG PT + tGoVi + tG V 4 G f + V 4 G f 

= tP + (l + tG )V 4 + tG Pf + (1 + tG )V 4 G f (31) 



with 



T = t 4 + T (32) 

This Eq.QJ]]) is identical to (0). 

Now the full breakup operator as given in (Q) is 

U = (1 + P)T + T 4 

= (1 + P)T + (1 + P)r 4 

= (1 + P)f (33) 

as required by (P~D|) . 

Finally the operator for elastic scattering QTOD can be rewritten, using ( P9"D and (2S) as 



£/■ = PG" 1 + PT + T 4 
= PG 1 + PT + r 4 

= PG 1 + PT + V 4 + V^GqT (34) 
which is identical to (P3"p. This concludes the verification of the equivalence of the old and 



new formulations. In practice we use the second form of Eq.fl3~3|) with the quantity r 4 to 
evaluate U, since it occurs naturally in Eq.(^) as an intermediate amplitude. 



We developed a computer code for the new forms (|21|), (|T^) and (|23|). This has significant 
practical advantages over the previous formulation. In the new approach one is simply 
iterating (|2T|), typically 18 times for total angular momentum and parity l/2 + followed by 
a Pade summation. The number of iterations for l/2~ is smaller and decreases further with 
increasing angular momentum. In the old scheme the set (|3|) expanded in power of V4 (see 



iT| ) had to be iterated. Since typically for l/2 + 10 powers of V4 had to be kept and for 
each order of V4 one has to iterate the corresponding Faddeev equations typically 13 times, 
the number of iterations in the new approach is drastically reduced. However, in the new 
approach the kernel in Eq.(pl|) includes V4 dependent parts and requires therefore more CPU 
time for its evaluation. Altogether we found in a realistic calculation a reduction of about a 
factor 4 in the CPU time for the new approach. 

We demonstrate now the feasibility of the new formalism and the numerical accuracy 
achieved with that form by showing a few 3N scattering observables evaluated with the 
old and new code. These are fullfledged calculations based on the NN potential CD Bonn 
(np) jn| kept up to j max = 2 (j is the total 2-body angular momentum) and the Tucson- 
Melbourne Ti — Ti exchange 3NF |J, also kept up to j max = 2. The cut-off parameter in the 
3NF has been chosen as A = 5.8 m^. 

Fig. H shows two spin observables in elastic nd scattering at 10.3 MeV. We have chosen 
cases where the effect of the 3NF is large. (The elastic differential cross section shows no 
effect of the 3NF at all.) The two curves with the 3NF included and evaluated in the old 
and new scheme overlap perfectly. This is generally true for all elastic observables. 

For the breakup process we show again two examples in Fig. |2] with large 3NF effects. 
One is a breakup cross section under np QFS conditions and the other a deuteron analyzing 
power under the same conditions. For that very small observable A y one can see tiny 
differences between the predictions of the old and new scheme. This is the worst case we 
found by checking in the standard breakup configurations cross sections and all analyzing 
powers, spin correlation coefficients and vector spin transfer coefficients. 

Summarizing, we have introduced a new single equation for 3N scattering including a 
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3NF. It is a direct generalisation of what has already been used for the 3N bound state. 
It is significantly faster to solve than our previous form used up to now |I2|| . We 
demonstrated the feasibility and accuracy of the new form. 
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FIGURES 

FIG. 1. Fig. 1: Spin correlation coefficient C yy and neutron to deuteron spin transfer coef- 
ficient if* v at E l £ h = 10.3 MeV. The solid curve is without 3NF and the long and short dashed 
curves are with 3NF evaluated with the old and new approach, respectively. 

FIG. 2. Fig. 2: The differential cross section and the deuteron analyzing power at E l $ b = 10.3 
MeV under np QFS conditions as a function of the scattering angle of one nucleon. Curves as in 
Fig. 1 
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